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ABSTRACT 
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q | We study the statistical isotropy (SI) of temperature fluctuations of the CMB as dis- 

tinct from Gaussianity. We present a detailed formalism of the bipolar power spectrum 
(T) ■ (BiPS) which was introduced as a fast method of measuring the statistical isotropy by 

Hajian & Souradeep 2003. The method exploits the existence of patterns in the real 
space correlations of the CMB temperature field. We discuss the applications of BiPS in 
constraining the topology of the universe and other theoretical scenarios of SI violation. 
Unlike the traditional methods of search for cosmic topology, this method is compu- 
tationally fast. We also show that BiPS is potentially a good tool to detect the effect 
| of observational artifacts in a CMB map such as non-circular beam, anisotropic noise , 

etc. Our method has been successfully applied to the Wilkinson Microwave Anisotropy 
,-£3 ■ Probe sky maps by Hajian et al. 2004, but no strong evidence of SI violation was found. 
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6 

52 ■ 1. Introduction 

03 

The first-year Wilkinson Microwave Anisotropy Probe ( WMAP) observations are consistent 
with predictions of the concordance ACDM model with scale-invariant and adiabatic fluctuations 
which have been generated during the inflationary epoch (Hinshaw et al. 2003; Kogut et al. 2003; 
Spergel et al. 2003; Page et al. 2003; 1). After the first year of WMAP data, the SI of the CMB 
anisotropy (i.e. rotational invariance of n-point correlations) has attracted considerable atten- 
tion. Tantalizing evidence of SI breakdown (albeit, in very different guises) has mounted in the 
WMAP first year sky maps, using a variety of different statistics. It was pointed out that the 
suppression of power in the quadrupole and octopole are aligned (Tegmark et al. 2004). Fur- 
ther "multipole-vector" directions associated with these multipoles (and some other low multi- 
poles as well) appear to be anomalously correlated (Copi et al. 2004; Schwarz et al. 2004). There 
are indications of asymmetry in the power spectrum at low multipoles in opposite hemispheres 
(Eriksen et al. 2004a; Hansen et al. 2004; Naselsky et al. 2004). Possibly related, are the results 
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of tests of Gaussianity that show asymmetry in the amplitude of the measured genus amplitude 
(at about 2 to 3c significance) between the north and south galactic hemispheres (Park 2004; 
Eriksen et al. 2004b; Eriksen et al. 2004c). Analysis of the distribution of extrema in WMAP sky 
maps has indicated non-gaussianity, and to some extent, violation of SI (Larson & Wandelt 2004). 
However, what is missing is a common, well defined, mathematical language to quantify SI (as 
distinct from non Gaussianity) and the ability to ascribe statistical significance to the anomalies 
unambiguously. 

In order to extract more information from the rich source of information provided by present 
(and future) CMB maps, it is important to design as many independent statistical methods as 
possible to study deviations from standard statistics, such as statistical isotropy (SI). Since SI can 
be violated in many different ways, various statistical methods can come to different conclusions. 
Because each method by design is more sensitive to a special kind of SI violation. 

As a statistical tool of searching for departures from SI, we use the bipolar power spectrum 
(BiPS) which is sensitive to structures and patterns in the underlying total two-point correlation 
function (Hajian & Souradeep 2003b; Souradeep Sz Hajian 2003). The BiPS is particularly sensi- 
tive to real space correlation patterns (preferred directions, etc.) on characteristic angular scales. 
In harmonic space, the BiPS at multipole £ sums power in off-diagonal elements of the covariance 
matrix, (ai m ai' m r) , in the same way that the 'angular momentum' addition of states Im, I'm' have 
non-zero overlap with a state with angular momentum \l — l'\ < I < I + 1'. Signatures, like a\ m and 
ai+nm being correlated over a significant range I are ideal targets for BiPS. These are typical of SI 
violation due to cosmic topology and the predicted BiPS in these models have a strong spectral 
signature in the bipolar multipole £ space (Hajian & Souradeep 2003a). The orientation indepen- 
dence of BiPS is an advantage since one can obtain constraints on cosmic topology that do not 
depend on the unknown specific orientation of the pattern (e.g., preferred directions). 

In the rest of this section we briefly introduce some basic statistical properties of CMB and 
discuss the statistical isotropy as distinct from Gaussianity. In sections 2 and 3 we introduce BiPS 
in real as well as harmonic space. In section 4 and 5 we define an unbiased estimator for BiPS 
and compute its cosmic variance in details. In section 6 we briefly introduce some possible sources 
of breakdown of SI. In section 7 we discuss the results of applying the method to the simulations 
and WMAP data. We draw our conclusions in section 8. For completeness we have also included 
details of calculations in 6 appendices to keep the paper easily readable. 

1.1. Statistical Properties of the CMB Anisotropy 

The temperature anisotropies of the CMB are described by a random field, AT(n) = T(fi) — To, 
on a 2-dimensional surface of a sphere (the so called surface of last scattering), where n = (6,4>) 
is a unit vector on the sphere and Tq = f ^ L T(h) represents the mean temperature of the CMB. 
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$ o(x) = / 7— 3$ (k)expik-x (2) 



The statistical properties of this field can be characterized by the n-point correlation functions 

(AT(ni)AT(n 2 ) • • • AT(n n )). (1) 

Here the bracket denotes the ensemble average, i.e. an average over all possible configurations of 
the field. If the field is Gaussian in nature, the connected part of the n-point functions disappears 
for n > 2. The non-zero (even-n)-point correlation functions can be expressed in terms of the 
2-point correlation function. As a result, a Gaussian distribution is completely described by the 
two-point correlation function. 

The statistics of the temperature fluctuations is Gaussian because they are linearly connected 
to quantum vacuum fluctuations of a very weakly interacting scalar field (the inflaton). Pertur- 
bations in primordial gravitational potential, $o> generate the CMB anisotropy, AT. The linear 
perturbation theory gives a linear relation between $o and AT. <3?o( x ) is assumed to be Gaus- 
sian, because in inflation, quantum fluctuations produce Gaussian random phase variations in the 
gravitational potential 

d 3 k 

(2^ 

with (<I>o(k)<I>Q(k / )) = P(k)5(k — k'), the power spectrum of potential fluctuations at that epoch. 
In the linear temperature approximation, in a spatially flat universe, the temperature anisotropies 
at the surface of last scattering can be linearly related to the potential fluctuations at that epoch 1 

/d 3 k 
__^-nr recMk)gT{kh (3) 

where r rec is the comoving conformal lookback time of the recombination epoch (r is the conformal 
time t = j dt/a), 4>(r) is a normalization constant for potential perturbations which is equal 
to unity today, and gr{k) is the linear radiation transfer function which describes the relationship 
between potential fluctuations and temperature fluctuations at recombination. For small k, gT{k) = 
1/3, which is known as the fact that for temperature fluctuations on super-horizon scales at the 
decoupling epoch, the Sachs- Wolfe effect (2) dominates and hence on large angular scales 

AT 1 

— (h) = -<S>Ls(nT rec ,T ) (4) 

On sub-horizon scales, gx oscillates (acoustic oscillation), and we need to evolve the coupled matter- 
radiation fluid, i.e. solve the Boltzmann photon transport equations coupled with the Einstein 
equations for gx- Codes like CMBFAST 2 compute gr(k) for different cosmological models. 

If we expand the gravitational potential into a series in powers of a density enhancement 

$(x,t) = qW + <S>W + ... (5) 



1 This description is very simple but naive. It does not take the line of sight effects such as integrated Sachs-Wolfe 
effect into account. For a more careful description see Appendix A. 

2 Available at http://www.cmbfast.org/ 
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The linear term $ 1 = <&o( x ) and is assumed to be Gaussian .Non-linearity in inflation makes $ 
weakly non-Gaussian and will lead to non-linear terms, such as <&( 2 \ in eqn. (5). However, since 
potential fluctuations in the early universe are small, ($ 2 ) ~ 1CT 9 , these non-linear effects are 
usually viewed as undetectable (Munshi et al. 1995; Spergel et al. 1999) (for a nice review on non- 
Gaussianity see (Bartolo et al. 2004)). Hence, here we limit our study to Gaussian CMB anisotropy 
field, where the two-point correlation function 

C(n, n ) = (AT(n)AT(n')) (6) 

contains all the statistical information encoded in the field. 

It is convenient to expand the temperature anisotropy field into spherical harmonics, the 
orthonormal basis on the sphere, as 

AT(n) = ^a im Y im (n) , (7) 

l,m 

where the complex quantities, a\ m are given by 

a lm = J dfi^(«)AT(n). (8) 

For a Gaussian CMB anisotropy, a\ m are Gaussian random variables and hence, the covariance 
matrix, {ai m a*, m ,}, fully describes the whole field. 



1.2. Statistical Isotropy of Temperature Anisotropies 

A random field 3>(r) is statistically homogeneous (in the wide sense) if the mean and the second 
moment (covariance) of that remain invariant under the coordinate transformation r — > r + Sr. 
That is to say, 

($(r)> = ($(r + Sr)), (9) 
C$(ri, r 2 ) = C$(ri + Sr, r 2 + Sr). 

The first condition implies that ($) = const. Putting Sr = — r 2 in the second condition we 
find that C$(ri, r 2 ) = C$(ri — r 2 , 0), i.e. the covariance of a statistically homogeneous field is 
dependent only on the difference r\ — r 2 , and not on ri and r 2 separately. So instead of C$(ri — r 2 , 0) 
we can simply write C$(r) which means: 

C*(r) = <*(n)$(n+r)>. (10) 

where $ = $ — ($). 

Statistically homogeneous fields, with C$(r) dependent only on the magnitude(but not on the 
direction) of the vector r = r 2 — ri connecting the points ri and r 2 

C*(r) = C*(r), r = |r| = s/F^E, (11) 
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are said to be statistically isotropic. The covariance of isotropic random fields is even and always 
real. 

For temperature anisotropies of the CMB, statistical isotropy is simply equivalent to rotational 
invariance of n-point correlation functions 

(AT(ftni)AT(ftn 2 ) • • • AT{TZh n )) = (AT(ni)AT(n 2 ) • • • AT(n n )), (12) 

where AT(JZn\) is the temperature anisotropy at lZh\ which is pixel h\ rotated by the rotation 
matrix TZ{a,(5,^f) for the Euler angles a,f3 and 7. For two point correlations of CMB, where 
l r i| = KiTred = const., the condition of statistical isotropy (eqn. 11) will translate into 

C(n, h') = C{6), 6 = arccos (n • n') (13) 

which implies that the two point correlation is just a function of separation angle, 6, between the 

two points on the sky. It is then convenient to expand it in terms of Legendre polynomials, 

^ 00 

C W = ^E( 2/ + 1 ) C '' P '( cos ^' ( 14 ) 
1=2 

where C\ is the widely used angular power spectrum. The summation over / starts from 2 because 
the I = term is the monopole which in the case of statistical isotropy the monopole is constant, 
and it can be subtracted out. The dipole I = 1 is due to the local motion of the observer and is 
subtracted out as well. 

To obtain the statistical isotropy condition in harmonic space, we expand the left hand side of 
eqn. (12) in terms of spherical harmonics, each AT(JZhi) may be expanded as 

AT(ftni) = J2 a hn Y lm(^i) (15) 

l,m 

l,m m' 

in which aj^ are harmonic coefficients in rotated coordinates and D l m , m (lZ) is the finite rotation 
matrix in the Zm-representation, Wigner D-function (Varshalovich et al. 1988). By substituting 
this into eqn. (12) and using the orthonormality law of spherical harmonics, eqn.(F2), we obtain the 
statistical isotropy condition in harmonic space 

E • • • a tmO D t im[ (K)Dt 2m , in)--. D 1 -^ (K) = (a hmi a hm2 • • • a lnm J. (16) 

m' 1 ,m,2---,m' n 

For a Gaussian field, only the covariance, {ai m a* llm ,), is important. As it is shown in Appendix B, 
the condition for statistical isotropy translates to a diagonal covariance matrix, 

( a lm a l' m ') = Ci5u'S mm i, (17) 

which means that the angular power spectrum of CMB anisotropy, Ci, tells us everything we need 
to know about the (Gaussian) CMB anisotropy. When (ai m a*, m ,) is not diagonal, the angular power 
spectrum C\ does not have all the information of the field and one should also take the off-diagonal 
elements into account. 
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1.3. Statistical Isotropy and Gaussianity 



Statistical isotropy and Gaussianity of CMB anisotropies are two independent and distinct 
concepts. The Gaussianity of the primordial fluctuations is a key assumption of modern cosmology, 
motivated by simple models of inflation. And since the statistical properties of the primordial 
fluctuations are closely related to those of the cosmic microwave background (CMB) radiation 
anisotropy, a measurement of non-Gaussianity of the CMB will be a direct test of the inflation 
paradigm. If CMB anisotropy is Gaussian, then the two point correlation fully specifies the sta- 
tistical properties. In the case of non-Gaussian fields, one must take the higher moments of the 
field into account in order to fully describe the whole field. Statistical isotropy is a condition which 
guarantees that statistical properties of CMB sky (e.g. n-point correlations) are invariant under 
rotation. Hence, in a SI model, the iso-contours of two point correlations where one point is fixed 
and the other scans the whole sky, are concentric circles around the given point. Any breakdown 
of statistical isotropy will make these iso-contours non-circular. In addition to cosmological mech- 
anisms, instrumental and environmental effects in observations easily produce spurious breakdown 
of statistical isotropy. 

A Gaussian field may or may not respect the statistical isotropy. In other words we can have 

1. Gaussian and SI models: In these models two point correlation and C\ have all the information 
of the field and C{n\,h-i) = C(h\ ■ 77,2) 

2. Gaussian but not SI models: two point correlation contains all the information of the field, 
but C(hi,n,2) / C(ni • r^)- So C\ is not adequate to describe the field and off-diagonal 
elements of covariance matrix should be taken into account. 

3. Non-Gaussian but SI models: two point correlation alone is not sufficient to describe the 
whole field and we need to take higher moments into account but every n-point correlation 
of the field is invariant under rotation. 

4. Non-Gaussian and non-SI models: neither two point correlation has all the information nor 
it is invariant under rotation. 

Here we confine ourselves to Gaussian fields where we only need to consider two point correlations. 
We will present a general way of describing fields of the kind (1) and (2) in the above list. 



Two point correlation of CMB anisotropies, C(h\, n.2), is a two point function on S 2 x S 2 , and 
hence can be expanded as 



2. The Bipolar Power Spectrum (BiPS) 




(18) 



h,l 2 ,L,M 
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where Afff are coefficients of the expansion (here after BipoSH coefficients) and {YJ 1 (hi)®Yi 2 {h 2 )}tM 
are the bipolar spherical harmonics which transform as a spherical harmonic with £, M with respect 
to rotations (Varshalovich et al. 1988) given by 

{Y ll (n 1 )®Y h (n 2 )} eM = £ CfZ lhm2 Y hmi (h 2 )Y hm2 (h 2 ), (19) 

raim2 

in which Cf^ il2m2 are Clebsch-Gordan coefficients. We can inverse-transform C(hi, n 2 ) to get the 
Aj^ 2 by multiplying both sides of eqn.(18) by {^' 1 (^i)<8)^^(«2)}|/m' an d integrating over all angles, 
then the orthonormality of bipolar harmonics, eqn. (F2), implies that 

AgE = J dn hl J dn A2 C{n u n 2 ) {Y h (ni) ® Y l2 (n 2 )}* m . (20) 

The above expression and the fact that C(ni, 7x2) is symmetric under the exchange of hi and h 2 
lead to the following symmetries of Af^f 

= (-l) ih+h - L) A^ 2 , (21) 
Ai M = Jiff 61^, A; = 0,1,2,-... 

We show in Appendix C that the Bipolar Spherical Harmonic (BipoSH) coefficients, Aj^f , are 
in fact linear combinations of off-diagonal elements of the covariance matrix, 

A-llh = X/ ( a 'i m i a «2rn 2 )( _ l) m2 Q £ imi/2-m 2 . ( 22 ) 



This means that Af^f completely represent the information of the covariance matrix. Fig. 1 shows 
how Af*f and Afff combine the elements of the covariance matrix. When SI holds, the covariance 
matrix is diagonal and hence 

A\f = (-lyc^l + l^SwSeoSMo, (23) 
A™ = (-l^v^T+lC,,^. 

BipoSH expansion is the most general way of studying two point correlation functions of CMB 
anisotropy. The well known angular power spectrum, C\ is in fact a subset of BipoSH coefficients. 

It is impossible to measure all Aj^f individually because of cosmic variance. Combining BipoSH 
coefficients into Bipolar Power Spectrum reduces the cosmic variance 3 . BiPS of CMB anisotropy 
is defined as a convenient contraction of the BipoSH coefficients 

«t= £|A^| 2 >0. (24) 

1,1', M 



3 This is similar to combining ai m to construct the angular power spectrum, Ci = J2 m \ a im\ 2 , to reduce the 
cosmic variance 
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The BiPS has interesting properties. It is orientation independent and is invariant under rotations 
of the sky. For models in which statistical isotropy is valid, BipoSH coefficients are given by 
eqn (23). And results in a null BiPS, i.e. Kg = for every positive £, 

ne = K 5 eo . (25) 

So, non-zero components of BiPS imply the break down of statistical isotropy. And this 
introduces BiPS as a measure of statistical isotropy. 

STATISTICAL ISOTROPY => k £ = W ^ 0. (26) 



3. BiPS Representation in Real Space 

Two point correlation of the CMB anisotropy is given by ensemble average over many universes 
(eqn. (6)). But in reality, there is only one universe and the ensemble average is meaningless unless 
the CMB sky is SI, where the two point correlation function C{6) can be well estimated by the 
average product of temperature fluctuations in two directions h and h' whose angular separation 
is 8 (see Fig. 2), this can be written as 

C{Q) = ( ^5(n ■ h' - cos 9) AT(h)AT(h') (27) 

J 4lT J 47T 

If the statistical isotropy is violated the estimate of the correlation function from a sky map 
given by a single temperature product 

C(ni,n 2 ) = AT(ni)AT(n 2 ) (28) 

is poorly determined. 

Although it is not possible to estimate each element of the full correlation function C(ni, n 2 ), 
some measures of statistical isotropy of the CMB map can be estimated through suitably weighted 
angular averages of AT(ni)AT(n 2 ). The angular averaging procedure should be such that the mea- 
sure involves averaging over sufficient number of independent measurements , but should ensure that 
the averaging does not erase all the signature of statistical anisotropy. Another important desirable 
property is that measure be independent of the overall orientation of the sky. Based on these con- 
siderations, we have proposed a set of measures of statistical isotropy (Hajian & Souradeep 2003b) 
which can be shown that is equivalent to the one in eqn. (24) 

k 1 = (2i + 1) 2 J dn m J dn n2 [ J dn x e (n) c(nh u nn 2 )} 2 . (29) 

In the above expression, C(TZh\, 7£n 2 ) is the two point correlation at 1Zh\ and 7£n 2 which are the 
coordinates of the two pixels hi and n 2 after rotating the coordinate system through an angle 00 
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where (0 < uj < it) about the axis n(0, $). The direction of this rotation axis n is defined by the 
polar angles G where (0 < < 7r) and where (0 < $ < 2-ir). x e is the trace of the finite rotation 
matrix in the £M-representation 

i 

x\n) = £ D e MM (n), (30) 

M=-i 

which is called the characteristic function, or the character of the irreducible representation of rank 
i. It is invariant under rotations of the coordinate systems. Explicit forms of X'iTV) are simple 
when 1Z is specified by a;, 0, <3?, then x^C^O is completely determined by the rotation angle uj and 
it is independent of the rotation axis n(0, $), 

X e (K) = xV) (31) 
sin [(21 + l)u;/2] 
sin [uj/2] 

And finally dIZ in eqn(29) is the volume element of the three-dimensional rotation group and is 
given by 

dIZ = 4 sin 2 | alto sin dO d<S> . (32) 
As it is shown in Appendix (D), this expression for BiPS may be simplified further as 

= J dn * J d&hCiquh) J dRxeiRMRq^Rfc) (33) 

For a statistically isotropic model C(ni, fi2) is invariant under rotation, and therefor C(1Zni, Ttfi'i) 
C(hi, fi2)- If we use this property to substitute in eqn.(33) for C(7Zhi, 7Zfi2), and if we use the 
orthonormality of x i { UJ ) ( see ec l n - we will recover the condition for SI, 

K e = K°5 eo . (34) 

The proof of equivalence of BiPS definition in eqns.(29) and (24) is as follows. In eqn.(29), we 
substitute (30) for x l {^) an d from 

IM' 

C(Tlh u TZh 2 ) = Y. A iXY. D MM'(mYh(ni) ® Y h (h 2 )} m , (35) 

hh M 

we substitute for C(R.ni, lZn.2) and use the orthonormality of D e MM (lZ) (see eqn.(F8)), we will 
obtain the eq.(24). Real-space representation of BiPS is very suitable for analytical compution of 
BiPS for theoretical models where we know the analytical expression for the two point correlation 
of the model, such as theoretical models in (Hajian & Souradeep 2003a). This formalism would 
immensely simplify analytical calculations. On the other hand, the harmonic representation of 
BiPS allows computationally rapid methods for BiPS estimation from a given CMB map. 
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4. Unbiased Estimator of BiPS 

In statistics, an estimator is a function of the known data that is used to estimate an unknown 
parameter; an estimate is the result from the actual application of the function to a particular set 
of data. Many different estimators may be possible for any given parameter. The above theory can 
be use to construct an estimator for measuring BipoSH coefficients from a given CMB map 

A™ = Y, Vwm>a lm a Vm , CZm' , (36) 

mm' 

where Wi is the Legendre transform of the window function. The above estimator is un-biased. 
Bias is the mismatch between ensemble average of the estimator and the true value. Bias for the 
BiPS is defined as ^Bg = {kg) — Kg and is equal to 

imi a hm' 1 ){ a hm 2 a hm' 2 ) + {a*l imi ai 2 m' 2 ){a*l 2 m2 a hm' 1 ) 

h,h mi,m' 1 m2,m' 2 

x / ^hmitem^hm'jhmi, ■ \°') 
M 

Therefore an un-biased estimator of BiPS is given by 

K £ = ^|i™| 2 -^, (38) 

ll'M 

The above expression for 25^ is obtained by assuming Gaussian statistics of the temperature fluc- 
tuations. The procedure is very similar to computing cosmic variance (which is discussed in the 
next section), but much simpler. However, we can not measure the ensemble average in the above 
expression and as a result, elements of the covariance matrix (obtained from a single map) are 
poorly determined due to the cosmic variance. The best we can do is to compute the bias for the 
SI component of a map 

_ t+h 

K^{Rf) sl = (2l + l)Y E C h C h W h W l2 [l + (-lYS hh ]. (39) 

h h=\e-h\ 

Note , the estimator R,£ is unbiased, only for SI correlation,i.e., (ki) = 0. Consequently, for 
SI correlation, the measured kg will be consistent with zero within the error bars given by a SI 
(Hajian & Souradeep 2003b). We simulated 1000 SI CMB maps and computed BiPS for them 
using different filters. The average BiPS of SI maps is an estimation of the bias which can be 
compared to our analytical estimation. Fig. 3 shows that the theoretical bias (computed from 
average C;) match the numerical estimations of average Kg of the 1000 realizations of the SI maps. 

It is important to note that bias cannot be correctly subtracted for non-SI maps. Non-zero 
ki estimated from a non-SI map will have contribution from the non-SI terms in full bias given in 
eq. (37). It is not inconceivable that for strong SI violation, 25^ over-corrects for the bias leading to 
negative values of kg. What is important is whether measured kg differs from zero at a statistically 
significant level. 
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Cosmic Variance 



Cosmic variance is defined as the variance of the estimator of the BiPS 

" 2 - 2 



a' 



< *t (40) 

We can analytically compute the variance of k f ~ using the Gaussianity of AT. Looking back at the 
eq.(29) we can see, we will have to calculate the eighth moment of the field 

(AT(n 1 )AT(n 2 )AT(n 3 )AT(n4)AT(n 5 )AT(n 6 )AT(n 7 )AT(n 8 )). (41) 

Assuming Gaussianity of the field we can rewrite the eight point correlation in terms of two point 
correlations. One can write a simple code to do that 4 . This will give us (8 — 1)!! = 7x5x3 = 105 
terms. These 105 terms consist of terms like: 

(AT(n 1 )AT(n 2 ))(AT(n 3 )AT(n4))(AT(n 5 )AT(n 6 ))(AT(n 7 )AT(n 8 )), (42) 

and all other permutations of them. On the other hand (Rf) has a 4 point correlation in it which 
can also be expanded versus two point correlation functions. If we form (Rj) — (Re) 2 , only 96 terms 
will be left which are in the following form 
,2£ + l 2 



8vr 2 



dfii---dfi 4 / dR / dR'xi{R)Xi{R')C{nuR'n 4 )C(n 2 ,R'7i 3 )C(Rh 1: n 4 )C(Rh 2 , n 3 ) 



(43) 

and all other permutations. In Appendix E we explain how we can handle these 96 terms to 
compute the following analytical expression for cosmic variance 

,(2^ + l) 2 



2Z + 1 



l:2l>£ 

e+h 

+ E E AC h CfM Wf 2 [{2i + 1) + Ff lh 
h h=\t-h\ 

2 



+ (-l)*(2* + l) + (l + 2(-l) 



(2£ + l) 2 



+ 16(-1) 



E 

l V .2h>e 



E c hW h 

h=\e-h\ 

e+h 

E, 
_, , , , 

h=\e-h\ 



C h C h W h W h- 



(44) 



Numerical computation of cr^ is fast. But the challenge is to compute Clebsch-Gordan coef- 
ficients for large quantum numbers. We use drc3j subroutine of netlib 5 in order to compute the 
Clebsch-Gordan coefficients in our codes. Again we can check the accuracy of our analytical estima- 
tion of cosmic variance by comparing it against the standard deviation of BiPS of 1000 simulations 
of SI CMB sky. The result is shown in Fig. 4 and shows a very good agreement between the two. 



4 F90 software implementing this is available from the authors upon request 
5 http: //www. netlib.org/slatec/src/ 
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6. Sources of Breakdown of Statistical Isotropy 

An observed map of CMB anisotropy, AT obs (n), is an n-dimensional vector (n = 12 x 512 2 for 
the WMAP data at HEALPix resolution N s id e = 512). This observed map contains the true CMB 
temperature fluctuations, AT(n), convolved with the beam and buried into noise and foreground 
contaminations. The observed map is related to the true map through this relation 

AT obs = BAT + n ( 45 ) 

in which B is a matrix that contains the information about the beam smoothing effect and n is the 
contribution from instrumental noise and foreground contamination. Hence, the observed map is a 
realization of a Gaussian process with covariance C = C T + C N + C res where C T is the theoretical 
covariance of the CMB temperature fluctuations, C N is the noise covariance matrix and C res is the 
covariance of residuals of foregrounds. Breakdown of statistical isotropy can occur in any of these 
parts. In general we can divide these effects into two kinds: 

• Theoretical effects: These effects are theoretically motivated and are intrinsic to the true 
CMB sky, AT. We discuss two examples of these effects, i.e. non-trivial cosmic topology and 
primordial magnetic fields, in the next sections. 

• Observational artifacts: In an ideally cleaned CMB map, the true CMB temperature fluctua- 
tions are completely extracted from the observed map. But this is not always true. Sometimes 
there are some artifacts (related to B or n) left in the cleaned map which may in principle 
violate the SI. These effects are explained in section 6.3. 



6.1. Signatures of Cosmic Topology 

The question of size and the shape of our universe are very old problems that have re- 
ceived increasing attention over the past few years (Ellis 1971; Lachieze-Rey k, Luminet 1995; 
Levin et al. 1998; Bond,Pogosyan & Souradeep 1998,2000; de Oliveira-Costa et al. 2003; de Oliveira-Costa et al. 
Dineen et al. 2004; Copi et al. 2004). Although a multiply connected universe sounds non-trivial, 
but there are theoretical motivations (Linde 2004; Levin 2002) to favor a spatially compact uni- 
verse. One possibility to have a compact flat universe is the consideration of multiply connected 
(topologically nontrivial) spaces. The oldest way of searching for global structure of the uni- 
verse is by identifying ghost images of local galaxies and clusters or quasars at higher redshifts 
(Lachieze-Rey & Luminet 1995). This method can probe the topology of the universe only on 
scales substantially smaller than the apparent radius of the observable universe. Another method 
to search for the shape of the universe is through the effect on the power spectrum of cosmic density 
perturbation fields. In principle, this effect can be observed in the distribution of matter in the 
universe and the CMB anisotropy. Over the past few years, many independent methods have been 
proposed to search for evidence of a finite universe in CMB maps. These methods can be classified 
in three main groups. 
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• Using the angular power spectrum of CMB anisotropics to probe the topology of the Uni- 
verse. The angular power spectrum, however, is inadequate to characterize the peculiar 
form of the anisotropy manifest in small universes of this type. It has been shown that 
a nontrivial topology breaks the SI down and as a result, there is more information in a 
map of temperature fluctuations than just the angular power spectrum (Levin et al. 1998; 
Bond,Pogosyan Sz Souradeep 1998,2000; Hajian Sz Souradeep 2003a). 

• The second class of methods are direct methods which rely on multiple imaging of the 
CMB sky. The most well known methods amongst these methods are S-map statistics 
(de Oliveira-Costa et al. 1996; de Oliveira-Costa et al. 2003) and the search for circles-in- 
the-sky (Cornish, Spergel Sz Starkman 1998). 

• Third class of methods are indirect probes which deal with the correlation patterns of the 
CMB anisotropy field by using an appropriate combination of coefficients of the harmonic 
expansion of the field (Dineen et al. 2004; Donoghue et al. 2004; Hajian Sz Souradeep 2003a; 
Copi et al. 2004). 

The BiPS method is one of the strategies for detecting observable signatures of a small compact 
universe. The BiPS method is essentially based on the search for correlation patterns in CMB. 
Using the fact that statistical isotropy is violated in compact spaces one could use the bipolar 
power spectrum as a probe to detect the topology of the universe. A simple example of application 
of BiPS in constraining topology of the universe is for a T 3 universe, where the correlation function 
is given by 

C(q, q') = L~ 3 Mkn) e^ 7 ^""^" 1 '*'), (46) 

n 

in which, n is 3-tuple of integers (in order to avoid confusion, we use q to represent the direction 
instead of h), the small parameter < 1 is the physical distance to the SLS along q in units of 
L/2 (more generally, L/2 where L = (L1L2-L3) 1 / 3 ) and L is the size of the Dirichlet domain (DD). 
When e is a small constant, the leading order terms in the correlation function eq. (46) can be 
readily obtained in power series expansion in powers of e. For the lowest wavenumbers |n| 2 = 1 in 
a cuboid torus 

C(q,q') « 2j]P*(27r/L i )cos(7reAAft) (47) 

i 

3 

« C l-e 2 |Ag| 2 + 3e 4 ^(A gi ) 4 

i=i 

where Aqi are the components of A.q = q — q' along the three axes of the torus and fa = L/Li- 
From this, the non-zero K£ can be analytically computed to be 

k 2 2 368 4 288 6 20736 8 

eg = 7r(1 " 4e + TT ~~ e + ^T e) 

«4 122887T 2 8 

cl ~ ~^f^ e (48) 
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K4 has the information of the relative size of the Dirichlet domain and one can use it to constrain 
the topology of the universe. A detailed description of determining the topology of the universe 
with BiPS is given in (Hajian & Souradeep 2003a). 

When CMB anisotropy is multiply imaged, the bipolar power spectrum corresponds to a 
correlation pattern of matched pairs of circles. Since it is orientation independent, the BiPS has to 
be computed only once for the CMB map irrespective of the orientation of the DD with respect to 
the sky. This method is very fast because there are fast spherical harmonic transform methods of 
the CMB sky maps. Using these methods it is very fast and efficient to compute the BiPS even for 
maps in WMAP resolution (HEALPix, N si d e = 512) in a few minutes on a single processor. We 
have used BiPS to put constraints on the topology of the Universe. The results are in preparation 
and will be reported soon(Hajian et al. 2004). 

The Aj^f signature, which is not discussed here, contains more details of orientation of the SI 
violation. It may be possible to detect the Aj^f too for strong violation of SI. 



6.2. Primordial Magnetic Fields 

It has been shown that a cosmological magnetic field, generated during an early epoch of 
inflation (Ratra 1992; Bamba et al. 2004), can generate CMB anisotropies (Durrer et al. 1998). 
The presence of a preferred direction due to a homogeneous magnetic field background leads to non- 
zero off-diagonal elements in the covariance matrix (Chen et al. 2004). This induces correlations 
between ai + \. m and a;_i jm multipole coefficients of the CMB temperature anisotropy field in the 
following manner 

{ai m a* Vm >) = $m,m'[5i,i>Ci + (<5; + i,/'-i + 5i-i,c+iA)] 1 (49) 

where D\ is the power spectrum of off-diagonal elements of the covariance matrix. For a Harrison- 
Peebles- Yu-Zel'dovich scale-invariant spectrum, D[ behaves as l~ 2 . More precisely, it is given by 

A = 4 x io- 16 r 2 (-^) 4 . (50) 

InG 

This clearly violates the statistical isotropy. If we substitute the above (a; m a z * m ,) into eqn. (22), 
we will get 

4t = W(-i)M2Ji + WA,iAo (51) 

+ Ai ^2 ( _1 ) miC, fimi«2-mi^i+i,«2-i 
ml=-ll,ll 

+ Ai ^2 ( -1 ) miC fimiJ2-mi^i-i,i2+i]- 
ml=-U,U 

The first line is the statistically isotropic part and it just contributes to kq. But the interesting 
parts are the two other lines. We can substitute this expression into eqn. (24) and compute the 
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BiPS predictions for magnetic fields. This work is still under process and will be reported in the 
near future (Hajian et al. 2004b). 



6.3. Observational Artifacts 

Foregrounds and observational artifacts (such as non-circular beam, incomplete/non-uniform 
sky coverage and anisotropic noise) would also manifest themselves as violations of SI. 

• Anisotropic noise : The CMB temperature measured by an instrument is a linear sum of the 
cosmological signal as well as instrumental noise. The two point correlation function then 
has two parts, one part comes from the signal and the other one comes from the noise 

C(ni, n 2 ) = C s (fn, n 2 ) + C N (h u n 2 ). (52) 

Both signal and noise should be statistically isotropic to have a statistically isotropic CMB 
map. So even for a statistically isotropic signal, if the noise fails to be statistically isotropic 
the resultant map will turn out to be anisotropic. The noise matrix can fail to be statistically 
isotropic due to non-uniform coverage. Also if the noise is correlated between different pixels 
the noise matrix could be statistically anisotropic. A simple example of this is the diagonal 
(but anisotropic) noise given by the following correlation 

C N (n,n')=a 2 (n)5 m ,. (53) 

This noise clearly violates the SI and will result a non-zero BiPS given by 

i 

m=— t 

where f^ m are spherical harmonic transform of the noise, fi m = f dilf,Y^ m (h)a 2 (h). 

• The effect of non-circular beam : In practice when we deal with data, it is necessary to take 
into account the instrumental response. The instrumental response is nothing but the beams 
width and the form of the beam and can be taken into account by defining a beam profile 
function B(h, n'). Here n denotes the direction to the center of the beam and n' denotes the 
direction of the incoming photon. The temperature measured by the instrument is given by 

AT (n) = J AT(h')B(h, ti)da h , (55) 

Using this relation to calculate the correlation function C(hi, n 2 ) = (AT (hi) AT (712)} one 
would get 

C(ni, n 2 ) = J dQn' J dn fl/l (AT(h / )AT(h"))B(h 1 , h')B(h 2 , h") (56) 
dn h > J dn iV ,C(h f ,h")B(h 1 , h')B(h 2 , h"). 
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Only for a circular beam where B(h, h') = B{n ■ n'), the correlation function is statistically 
isotropic, C(h\, n 2 ) = C{n\ ■ fi2)- Breakdown of SI is obvious since even C\ get mixed 
for a non-circular beam, C\ = Ylv ^wCy (Mitra et al. 2004). Non-circularity of the beam 
in CMB anisotropy experiments is becoming increasingly important as experiments go for 
higher resolution measurements at higher sensitivity. 

• Mask effects : Many experiments map only a part of the sky. Even in the best case, contam- 
ination by galactic foreground residuals make parts of the sky unusable. The incomplete sky 
or mask effect is another source of breakdown of SI. But, this effect can be readily modeled 
out. The effect of a general mask on the temperature field is as follows 

AT masked (n) = AT(n)W(n), (57) 

where W(h) is the mask function. One can cut different parts of the sky by choosing appropri- 
ate mask functions. Masked a; m coefficients can be computed from the masked temperature 
field, 



^masked 
a lm 



= J &T masked (n)Y* m (n) d£ln (58) 
= ^ahmi J Yhmi(n)Y^(h)W(h)dnn. 

l\m\ 

Where a; imi are spherical harmonic transforms of the original temperature field. We can 
expand W(h) in spherical harmonics as well 

W{n) = J2 w imYim(n), (59) 

Im 

and after substituting this into eqn. (58) it is seen that the masked a\ m is given by the effect 
of a kernel K 1 ^ 1 on original a lm (Prunet et al. 2004) 



masked 
a lm 



E^C- (60) 

hmi 



The kernel contains the information of our mask function and is defined by 

Kim 1 = j Y hmi {h)Y hm MYL{n)dQ h (61) 




(2Z 1 + 1)(2Z 2 + 1) , lm 

Air(2l + 1) Ly li0l20 Ly hm 1 l2m 2 - 

In the last step of the above expression we used rule F3 of spherical harmonics to simplify 
the expression. The covariance matrix of a masked sky will no longer have the diagonal form 
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of eqn.(17) because of the action of the kernel 



/ masked masked * 
\ a lm a l'm' 



) 




-hrm j^l'i m 'i 
Im I'm' 



(62) 



-hmi Tf l 'i m i 
Im "'I'm' 



-hmi 
I'm' ■ 



li,mi 



This clearly violates the SI and results a non-zero BiPS for masked CMB skies. In the next 
section we apply a galactic mask to ILC map and show that signature of this mask on BiPS 
is a rising tail at low £, {I < 20). 

• Residuals from foreground removal : Besides the cosmological signal and instrumental noise, a 
CMB map also contains foreground emission such as galactic emission, etc. The foreground is 
usually modeled out using spectral information. However, residuals from foreground subtrac- 
tions in the CMB map will violate SI. Interestingly, BiPS does sense the difference between 
maps with grossly different emphasis on the galactic foreground. As an example we compute 
BiPS of a Wiener filtered map in the next section to make this point clear (see Fig. 8). It 
shows a signal very similar to that of a galactic cut sky. This can be understood if one writes 
the effect of the Wiener filter as a weight on some regions of the map 



This explains the similarity between a cut sky and a Wiener filtered map. The effect of 
foregrounds on BiPS still needs to be studied more carefully. 



We carry out measurement of the BiPS, on the following CMB anisotropy maps 

A) a foreground cleaned map (denoted as 'TOH') (Tegmark et al. 2004), 

B) the Internal Linear Combination map (denoted as 'ILC in the figures) (Bennett et al. 2003), 
and 

C) a customized linear combination of the QVW maps of WMAP with a galactic cut (denoted 



AT W 



(n) = AT(h)(l + W(n)). 



(63) 



7. 



BiPS Results from WMAP 



as 'CSSK'). 



Also for comparison, we measure the BiPS of 



D) a Wiener filtered map of WMAP data (denoted as Wiener') (Tegmark et al. 2004), and 

E) the ILC map with a 10° cut around the equator (denoted as 'Gal. cut.'). 
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Angular power spectra of these maps are shown in Fig. 5. The best fit theoretical power spectrum 
from the WMAP analysis 6 (Spergel et al. 2003) is plotted on the same figure. C\ from observed 
maps are consistent with the theoretical curve, Cf , (except for low multipoles. The bias and 
cosmic variance of BiPS depend on the total SI angular power spectrum of the signal and noise 
C\ = Cf + . However, we have restricted our analysis to / < 60 where the errors in the 
WMAP power spectrum is dominated by cosmic variance. It is conceivable that the SI violation 
is limited to particular range of angular scales. Hence, multipole space windows that weigh down 
the contribution from the SI region of multipole space will enhance the signal relative to cosmic 
error, a SI . We use simple filter functions in I space to isolate different ranges of angular scales; a 
low pass, Gaussian filter 

W?(l s ) = exp(-(/ + l/2) 2 /(Z s + 1/2) 2 ) (64) 
that cuts off power on small angular scales (< l/l s ) and a band pass filter, 

Wf(l t , l s ) = [2(1 - J ((/ + l/2)/(l t + 1/2)))] exp(-(Z + l/2) 2 /(/ s + 1/2) 2 ) (65) 

that retains power within a range of multipoles set by k and l s . The windows are normalized 
such that + l/2)/(Z(Z + l))Wi = 1, i.e., unit rms for unit flat band power Q = 1/(1(1 + 1)). 
The window functions used in our work are plotted in figure 5. We use the Cf to generate 1000 
simulations of the SI CMB maps. a; m 's are generated up to an l ma x of 1024 (corresponding to 
HEALPix resolution N s ^ e = 512). These are then multiplied by the window functions Wp(l s ) and 
W l (It, l s ). We compute the BiPS for each realization. Fig.5 shows that the average power spectrum 
obtained from the simulation matches the theoretical power spectrum, Cf, used to generate the 
realizations. We use Cf to analytically compute bias and cosmic variance estimation for K£. This 
allows us to rapidly compute BiPS with la error bars for different theoretical Cj . 

We use the estimator given in eqn.(36) to measure BiPS for the given CMB maps. We compute 
the BiPS for all window functions shown in Fig 5. Results for three of these windows are plotted 
in Figs. 6, 7, and 8. In the low-/ regime, where we have kept the low multipoles, BiPS for all three 
given maps are consistent with zero (Fig. 6). But in the intermediate-/ regime (Fig. 7), although 
BiPS of ILC and TOH maps are well consistent with zero, the CSSK map shows a rising tail in BiPS 
due to the galactic mask. To confirm it, we compute the BiPS for the ILC map with a 10-degree 
cut around the galactic plane (filtered with the same window function) . The result is shown on the 
top panel of Fig. 7. Another interesting effect is seen when we apply a Wp (20,45) filter, where 
Wiener filtered map has a non zero BiPS very similar to that of CSSK but weaker (Fig. 8). The 
reason is that Wiener filter takes out more modes from regions with more foregrounds since these 
are inconsistent with the theoretical model. As a result, a Wiener filtered map at Wj s (20, 45) filter 
has a BiPS similar to a cut sky map. The fact that Wiener map has less power at the Galactic 
plane can even be seen by eye! Comparing Fig. 7 to Fig.8 we see that using different filters allows 



6 Based on an LCDM model with a scale-dependent (running) spectral index which best fits the dataset comprised 
of WMAP, CBI and ACBAR CMB data combined with 2dF and Ly-a data 
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us to uncover different types of violation of SI in a CMB map. In our analysis we have used a set 
of niters which enables us to probe SI breakdown on angular scales I < 60. 

The BiPS measured from 1000 simulated SI realizations of Cj is used to estimate the prob- 
ability distribution functions (PDF), p(k(). A sample of the PDF for two windows is shown in 
Fig. 9. Measured values of BiPS for ILC, TOH and CSSK maps are plotted on the same plot. BiPS 
for ILC and TOH maps are located very close to the peak of the PDF. We compute the individual 
probabilities of the map being SI for each of the measured k,£. This probability is obtained by 
integrating the PDF beyond the measured k£. To be precise, we compute 



The probabilities obtained are shown in Figs. 10 and 12 for W s (20, 30), W Z G (40) and W t G {4). 
The probabilities for the Wf(20, 30) window function are greater than 0.25 and the minimum 
probability at ~ 0.05 occurs at K4 for W G (40). The reason for systematically lower SI probabilities 
for Wi (20, 30) as compared to Wi G (40) is simply due to lower cosmic variance of the former. The 
contribution to the cosmic variance of BiPS is dominated by the low spherical harmonic multipoles. 
Filters that suppress the a\ m at low multipoles have a lower cosmic variance. 

It is important to note that the above probability is a conditional probability of measured ki 
being SI given the theoretical spectrum Cf (used to estimate the bias). A final probability emerges 
as the Bayesian chain product with the probability of the theoretical Cf used given data. Hence, 
small difference in these conditional probabilities for the two maps are perhaps not necessarily 
significant. Since the BiPS is close to zero, the computation of a probability marginalized over 
the Cf may be possible using Gaussian (or, improved) approximation to the PDF of K£. The 
important role played by the choice of the theoretical model for the BiPS measurement is shown 
for a W\ that retains power in the lowest multipoles, I = 2 and I = 3. Assuming Cf, there 
are hints of non-SI detections in the low €s (left panel of Fig. 11). We also compute the BiPS 
using a Cf for a model that accounts for suppressed quadrupole and octopole in the WMAP 
data (Shafieloo & Souradeep 2004). The mild detections of a non zero BiPS vanish for this case 
(right panel of Fig. 11). 



The SI of the CMB anisotropy has been under scrutiny after the release of the first year of 
WMAP data. We use the BiPS which is sensitive to structures and patterns in the underlying 
total two-point correlation function as a statistical tool of searching for departures from SI. We 
carry out a BiPS analysis of WMAP full sky maps. We find no strong evidence for SI violation 
in the WMAP CMB anisotropy maps considered here. We have verified that our null results 




(66) 




8. Discussion and Conclusion 
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are consistent with measurements on simulated SI maps. The BiPS measurement reported here 
is a Bayesian estimate of the conditional probability of SI (for each K£ of the BiPS) given an 
underlying theoretical spectrum Cf . We point out that the excess power in the Cf with respect 
to the measured C/ from WMAP at the lowest multipoles tends to indicate mild deviations from 
SI. BiPS measurements are shown to be consistent with SI assuming an alternate model Cj that 
is consistent with suppressed power on low multipoles. Note that it is possible to band together 
Kg measurements to tighten the error bars further. The full sky maps and the restriction to low 
/ < 60 (where instrumental noise is sub-dominant) permits the use of our analytical bias subtraction 
and error estimates. The excellent match with the results from numerical simulations is a strong 
verification of the numerical technique. This is an important check before using Monte-Carlo 
simulations in future work for computing BiPS from CMB anisotropy sky maps with a galactic 
mask and non uniform noise matrix. 

There are strong theoretical motivations for hunting for SI violation in the CMB anisotropy. 
The possibility of non-trivial cosmic topology is a theoretically well motivated possibility that has 
also been observationally targeted (Ellis 1971; Lachieze-Rey & Luminet 1995; Levin 2002; Linde 2004). 
The breakdown of statistical homogeneity and isotropy of cosmic perturbations is a generic fea- 
ture of ultra large scale structure of the cosmos, in particular, of non trivial cosmic topology 
(Bond,Pogosyan & Souradeep 1998,2000). The underlying correlation patterns in the CMB anisotropy 
in a multiply connected universe is related to the symmetry of the Dirichlet domain. The BiPS 
expected in flat, toroidal models of the universe has been computed and shown to be related 
to the principle directions in the Dirichlet domain (Hajian & Souradeep 2003a). As a tool for 
constraining cosmic topology, the BiPS has the advantage of being independent of the overall 
orientation of the Dirichlet domain with respect to the sky. Hence, the null result of BiPS can 
have important implication for cosmic topology. This approach complements direct search for sig- 
nature of cosmic topology (Cornish, Spergel & Starkman 1998; de Oliveira-Costa et al. 1996) and 
our results are consistent with the absence of the matched circles and the null S-map test of the 
WMAP CMB maps (Cornish et al. 2003; de Oliveira-Costa et al. 2003). Full Bayesian likelihood 
comparison to the data of specific cosmic topology models is another approach that has applied 
to COBE-DMR data (Bond,Pogosyan & Souradeep 1998,2000). Work is in progress to carry out 
similar analysis on the large angle WMAP data. We defer to future publication, detailed ana- 
lyzes and constraints on cosmic topology using null BiPS measurements, and the comparison to 
the results from complementary approaches. There are also other theoretical scenarios that pre- 
dict breakdown of SI that can be probed using BiPS, e.g., primordial cosmological magnetic fields 
(Durrer et al. 1998; Chen et al. 2004). 

The null BiPS results also has implications for the observation and data analysis techniques 
used to create the CMB anisotropy maps. Observational artifacts such as non-circular beam, inho- 
mogeneous noise correlation, residual stripping patterns, etc. are potential sources of SI breakdown. 
Our null BiPS results confirm that these artifacts do not significantly contribute to the maps studied 
here. Foreground residuals can also be sources of SI breakdown. The extent to which BiPS probes 
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foreground residuals is yet to be fully studied and explored. We do not see any significant effect of 
the residual foregrounds in ILC and the TOH maps as it was mentioned by (Eriksen et al. 2004c). 
This can not be necessarily called a discrepancy between the two results unless we know what 
should have been seen in the BiPS. The question is if the signal is strong enough and whether the 
effect smeared out in bipolar multipole space within our angular /-space window. On the other 
hand, the very fact that BiPS does show a strong signal for the Wiener filtered map, mean that at 
some level BiPS is sensitive to galactic residuals. 

In summary, we study the Bipolar power spectrum (BiPS) of CMB which is a promising 
measure of SI. We find null measurements of the BiPS for a selection of full sky CMB anisotropy 
maps based on the first year of WMAP data. Our results rule out radical violation of statistical 
isotropy, and are consistent with null results for matched circles and the S-map tests of SI violation. 
Pending a more careful comparison, the results do not necessarily conflict with a number of other 
statistical tests, that hint at SI violation at low to modest statistical significance. 

The authors are very thankful for close interaction with David Spergel, Glenn Starkman and 
Neil Cornish. We also acknowledge useful discussions with Dick Bond, Dmitry Pogosyan and Carlo 
Contaldi. AH thanks Istvan Szapudi, Simon Prunet, Uros Seljak and Jacques Delabrouille for 
illuminating discussions. Use of the HPC facility of IUCAA is acknowledged. 



A. More on Gaussianity 

The temperature anisotropy at a given space-time point (x, r) can be written as a superposition 
of plane waves in A;-space 

AT(x, n, t) = j d 3 k e ik x AT(k, n, r) . (Al) 

As discussed in (Ma & Bertschinger 1995), the dependence on n arises only through k-n. Therefore 
AT(k, n, r) may be represented as 

oo 

AT(k, n,r) = Y / (-i) l (2l + 1) AT, (k, t)P, (k ■ n) . (A2) 

1=0 

The anisotropy coefficients ATJ (k, r) are random variables whose amplitudes and phases depend 
on the initial perturbations and can be written as 

AT,(k,r)=^o(k)Ar,(fc,r), (A3) 

where 0o(k) is the initial potential perturbation and AT/(/c,r) is the photon transfer function, i.e. 
solution of Boltzmann equation with 0o(k) = 1. In general, ATj(fc,r) is given by the integral 
solution 

AT l ((3)= dr^(T -r)S(P,T), (A4) 
J o 
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in which S(j3, r) is the source function, is the ultra-spherical Bessel function and (3 = k — K, 
where K = Hq(Qo — 1) is the curvature (Zaldarriaga et al. 1998). In the flat case, the above 
expression has a very simple form 

AT(k)= H drj^-^Sfcr), (A5) 
Jo 

and hence, the temperature anisotropics are given by 

AT(x,n) = j d 3 ke ik *Mk)[j^ dTe ih >r (T - T0 *>S(k,T)}. (A6) 

Equation (A6) clearly shows that Gaussian perturbations in primordial gravitational potential 
result a Gaussian CMB temperature anisotropy field. 



B. SI condition in harmonic space 

We substitute the spherical harmonic coefficients, a; m , from eqn. (8) in the covariance matrix 
to express it in terms of the two point correlation 

(a lm a* Vm ,) = j j da^n>Y; m (m>nAn')(AT(h)AT(ti)) (Bl) 

= J j &n f MniYC m {m> m >(n')C(h,ti). 

Under the condition of SI, eqn. (13), and after making use of the spherical harminc expansion of 
Legendre polynomials, eqn. (F6), in eqn. (14), the above expression will become 

{ai m a* Vm ,) = ! I 'd^dn^^y^^O^Q^^^y^^n') (B2) 

J J h mi 

= y^Qi (5// 1 Si'h s mmi s m ' mi 

h m i 

= Ci5u'5 mm >, 

which is the SI condition in harmonic space, eqn. (17), and in the last line, the orthonormality 
relation of spherical harmonics, eqn. (F2) has been used. 



C. Mii 2 are l mear combinations of (ai m ai, m ,) 

The BiPoSH coefficients, Aj^, are linear combinations a\ m . To see that, in the definition of 
two point correlation we can substitute for AT/T from eqn. (7), which will lead to the expansion 
of C(hi,ri2) in terms of spherical harmonics 

C(ni,n 2 ) = <AT(ni)Ar(n2)) (CI) 

= Y( a lrriaVm>) Y lrri{ni)Y{; m \h2). 
Im I'm' 



This can be substituted in the definition of Affi and eqn. (20) can be written as 

4Y 2 = jdSlij d^Cih^hiWiAhi) x Y l2 (n 2 m M (C2) 

= IdSlJ dQ 2 Y J Y.^4rn')yirn{ni)Y l J ml {h^ £ C^^Yf^^Y^h,) 
Im I'm 1 m 1 m 2 

= IdSlJ dn 2 Y,Y,( a ^rn')Yi m (n l )Y i r n An2) E Q^ li2m2 ^ mi (ni)(-l) m2 H 2 - m2 (n 2 ) 

Im I'm' m\m,2 
= E ( _1 ) m2 ( a hm 1 0'* 2 - m2 }Ct 1 m 1 hr 



mim 2 

2^2 



By changing the dummy variable 771,2 — ► -m2 we will get 

■^ilh = X] ( a 'i"ii a (2m 2 )( _1 ) m2 ^mii 2 -m 2 - (C3) 



mim 2 



D. Simplifying the real-space expression for BiPS 

The BiPS in real-space is given by the following expression. 

"t = (^r) 2 / dn i / dn 2 f dR lX t(Ri) j dR 2 xe(R2)C(R 1 quR 1 q2)C(R2quR 2 q2) (Dl) 

Now if we change the variables in the following way 

R lQl = fa (D2) 

Rlfa = fa 

We will have 

*t = (^J^) 2 / ^1 / j dRiXt(Ri) J dR 2 xe(R2)C(q 1 ,fa)C(R 2 R^ 1 qi,R2Ri 1 fa), (D3) 
and with 

i? = i?^- 1 (D4) 

we will get 

= j dQl J d ^ j dR iXe(Ri) j dR Xe (RRi)C(fa,fa)C(Rfa,Rfa). (D5) 

Using eqn. (F9) we can write 

J dR lX e(Ri)xe(RRi) = EE f d Ri D ^rnARi) D L m '(R)D e mlm (Ri) (D6) 

m 1 mm f 
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mi mm' 

8vr 2 



^ ^ ] ^mimi (^) 



2£ + 

mi 

8 ' 2 *<*). 



2£ + 1 

So, the K£ can be written in this form: 

2^ + 1 



8vr 2 



y dfii J dfi 2 C(gi,52) y dR X t(R)C(Rqi,Rq2). (D7) 



E. More on Cosmic Variance 

The cosmic variance, (k 2 ) — (re^) 2 , has 96 terms similar to the one shown in eqn. (43) but 
with different combinations of R's and re's. Among them, 40 terms contain at least one term like 
C(Rhi, Rh 2 ). These terms only contribute to kq because when SI holds, C(Rhi, Rfi2) = C(h\,n 2 ) 
and hence the integral in eqn. (D7) will be proportional to 5iq, which is not in the interest of us. 
So, we can ignore them and we will be left with 56 terms. These terms, depending on the number 
of connected correlations in them, can be classified into 4 major groups: 

1. One- loop terms, such as 

C(rei, R'n 3 )C(n 3 , Rh 2 )C(h 2 , R'h A )C{h A , Rhi). (El) 

These terms get reduced to expressions with a Cf factor and 4 Clebsch-Gordan coefficients 
(summation symbols are omitted for brevity), 

There are 24 terms like this but with different combinations of indices. 

2. Two-loop terms - type I, such as 

C(ni, n 3 )C(R'h 3 , Rh 2 )C(h 2 , Rhi)C{n A , R'n A ). (E3) 
These terms get reduced to expressions with a Cf Ci 2 factor and 4 Clebsch-Gordan coefficients 

Zi '7 h—mih—ms li—mslimz l\m\l-i—mi l\—mzh—m-r' V 1 -^/ 

There are 16 terms like this but with different combinations of indices. 

3. Two-loop terms - type II, such as 

C(ni, R'n 4 )C(n 4 , Rhi)C{n 2 , R'n 3 )C(Rh 2 , n 3 ). (E5) 
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These terms get reduced to expressions with a Cf Cf factor and 4 Clebsch-Gordan coefficients 

h I3 h— mils— rnz^ Ixm^lzm-j^ l\m$lzrn~i m\lz — 1713 • K*- 1 ®) 

There are 8 terms like this but with different combinations of indices. 

4. Three-loop terms, such as 

C(ni, h 3 )C(R'h 3 , Rh 1 )C{h 2 , Rh 2 )C(h 4 , R'h 4 ). (E7) 

These terms get reduced to expressions with a Cf Ci 2 Ci 3 factor and 4 Clebsch-Gordan coeffi- 
cients 

n 2 n n n tM n lM n lM' n tM' (T7-o\ 

There are 8 terms like this but with different combinations of indices. 

Using symmetry properties of Clebsch-Gordan coefficients, eqn. (F10), and their summation 
rules, eqn. (Ffl), it is possible to simplify these 56 terms immensely. After tedious algebra we 
would obtain the final expression for the kernel of cosmic variance of the bipolar power spectrum: 



4 clCl\(2i + l)+F[ il2 ]{hl 2 l} + 

8 c ^ ChCl3 ^Ti' {hhe} {hhe} + 



16 {-ifcldj^^- {hh£}{hh£} (E9) 



2h + l 

where the 3j-symbol {l\h£} is given by 



{hhi} = 



1 if h + h + £ is integer and \h - h\ < i < (h + h), 
otherwise, 



and Ff il3 is defined as 



h h t 

pt \ ^ \ ^ \ ^ s-itM plM ptM' plM' fFlfll 

hh / j / j / j Ii—mil3—rn3 lim2l3rn4 l3rn4hrni l3—rn3h—m2' \I2j±v) 

m\m,2=—l\ 77137714=— £3 M,M'=—£ 

This should be summed over all indices except £ to result the analytical expression for the cosmic 
variance given in eqn. (44). 
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F. Useful Standard Integrals and Summation Rules 

For completeness we list a set of standard integrals and summation rules which are needed in 
BiPS calculations. Orthonormality of spherical harmonics 

J dn h Y em (h) = V4^5 eo Smo, (Fl) 

(F2) 

J dn h Y hm Am 2 rn 2 (n)K m3 (n) = ^ to^CW ( F3 ) 

(F4) 

Symmetry property of spherical harmonics 

YLin) = (-l) m Y^ m (h). (F5) 
Spherical harmonic expansion of Legendre polynomials 



47T 



Pl(n ■ n') = £ YS n (A)Y lm (A')- (F6) 

m=—l 

Orthonormality relation of x^( w ) 

£ X h (u) X h M sin 2 ^dco = |«5^ 2 . (F7) 

Orthonormality of D e MM (TZ) 

J d<RDf IMI (n)D l mml (n) = ^^6 u S mM 6 m/MI , (F8) 
]T D e mmf (Ri)Di, m (R 2 ) = X e(RiR2) (F9) 

mm' 

Symmetry properties of Clebsch-Gordan coefficients 

= (- 1 ) a+6 " c c^' ( F1 °) 

/-icy _ /_-. Na+b_ c ^,c-7 

Summation rules of Clebsch-Gordan coefficients 



J2 C albp C aab/3 = Scc'tjj' {abc}{abc'} 

a/3 



Y. C ZbfZbV = Striate} (Fll) 



-27- 



REFERENCES 

Bamba, K., Yokoyama, J., Phys.Rev. D69 (2004) 043507. 
Bardeen, J. M., Steinhardt, P. J. & Turner, M. S. 1983, Phys. Rev. D 28, 679. 
Bartolo, N., Komatsu, E., Matarrese, S., Riotto, A., preprint (astro-ph/0406398) 
Bennett, C. L., et.al., 2003, Astrophys. J. Suppl., 148, 1. 

Bond, J. R. et al. 2004, Int. J. Theor. Phys. 2004, ed. Verdaguer, E., "The Early Universe: Con- 
fronting theory with observations" (June 21-27, 2003) (astro-ph/0406195). 

Bond, J. R., Pogosyan, D. & Souradeep,T. 1998, Class. Quant. Grav. 15, 2671; ibid. 2000, Phys. 
Rev. D 62,043005;2000, Phys. Rev. D 62,043006. 

Chen, G., Mukherjee, P., Kahniashvili, T., Ratra, B., Wang, Y., Astrophys.J. 611 (2004) 655. 

Coles, P. et al. 2003, preprint (astro-ph/0310252). 

Copi, C. J., Huterer, D. & Starkman, G. D. 2004, Phys. Rev. D. in press, (astro-ph/0310511). 

Cornish, N.J., Spergel, D.N. & Starkman, G. D. 1998, Class. Quantum Grav., 15, 2657 

Cornish, N. J., Spergel, D., Starkman, G. , Komatsu, E., 2004, Phys.Rev.Lett. 92, 201302. 

Cruz, M., Martinez-Gonzalez, E., Vielva, P., Cayon, L., preprint (astro-ph/0405341). 

de Oliveira-Costa, A., Tegmark, M., Zaldarriaga, M. & Hamilton, A. 2004, Phys. Rev.D69, 063516. 

de Oliveira-Costa, A. Smoot, G. F., Starobinsky, A. A., 1996, ApJ 468, 457. 

Dineen, P. , Rocha, G., Coles, P, preprint (astro-ph/0404356). 

Donoghue, E. P., and Donoghue, J. F., preprint (astro-ph/0411237). 

Durrer, R., Kahniashvili, T. and Yates, A., 1998, Phys. Rev. D 58, 3004. 

Ellis, G. F. R. 1971, Gen. Rel. Grav. 2, 7. 

Eriksen, H. K. et al., 2004, Astrophys. J 605, 14. 

Eriksen, H. K. et al., 2004, Astrophys. J. 612, 64. 

Eriksen, H. K. et al., 2004, Astrophys. J. 612, 633. 

Gaztanaga, E. & Wagg, J. 2003, Phys. Rev. D68 021302. 

G6rski,K. M., Hivon, E., Wandelt, B. D. 1999, in "Evolution of Large-Scale Structure", eds. A.J. 
Banday, R.S. Sheth and L. Da Costa, PrintPartners Ipskamp, NL, pp. 37-42 (also astro- 
ph/9812350). 



-28- 



Guth, A. H. & Pi, S.-Y. 1982, Phys. Rev. Lett., 49, 1110. 
Hajian, A. & Souradeep, T., 2003a preprint (astro-ph/0301590). 
Hajian, A. and Souradeep, T., 2003b, ApJ 597, L5 (2003). 
Hajian, A. & Souradeep, T., 2005, ApJ 618, L63. 

Hajian, A., Pogosyan, D. Souradeep, T., Contaldi, C, Bond, R., 2004 in preparation; Proc. 20th 
IAP Colloquium on Cosmic Microwave Background physics and observation, 2004. 

Hajian, A., Chen, G., Souradeep, T., Kahniashvilli, T., Ratra, B., 2004 in preparation. 

Hansen, F. K., Banday, A. J. & Gorski, K. M. 2004, preprint, (astro-ph/0404206) 

G. Hinshaw, Astrophys. J. Suppl.,(2003) ,148, 135. 

Komatsu, E. et.al., 2003, ApJS, 148, 119. 

Kogut A. et al., 2003, Astrophys. J. Suppl.,148, 161. 

Lachieze-Rey, M. and Luminet,J. -P. 1995, Phys. Rep. 25, 136. 

Larson, D. L. & Wandelt, B. D. 2004, preprint, (astro-ph/0404037). 

Levin, J. 2002, Phys. Rep. 365, 251. 

Levin J., Scannapieco E., Silk J., (1998), Class. Quant. Grav. 15, 2689. 

Linde, A., (2004), JCAP 0410, 004, (astro-ph/0408164). 

Luminet, J.-P. et al. 2003, Nature 425 593. 

Ma, C.-R, & Bertschinger, E. 1995, ApJ, 455, 7 

Mitra, S., Sengupta, A., S., Souradeep, T., preprint astro-ph/0405406 

Munshi, D., Souradeep, T., and Starobinsky, A. A., (1995), Astrophys. J., 454, 552. 

Page L. et al, 2003, Astrophys. J. Suppl.,148, 233. 

Naselsky et al, 2004, preprint, (astro-ph/0405523); ibid (astro-ph/0405181), 2003; ApJ,599,L53 and 
references therein. 

Park, C, 2004, Mon. Not. Roy. Astron. Soc. 349, 313. 

Peiris H.V. et al., 2003, Astrophys. J. Suppl.,148, 213. 

Prunet, S., Uzan, J., Bernardeau, F., Brunier, T., 2004, preprint (astro-ph/0406364). 



-29- 



Ratra, B., (1992), Astrophys.J. 391, LI. 

Sachs, R. K., and Wolfe, A. M. 1967, ApJ, 147, 73. 

Schwarz, D. J. et al., 2004, preprint (astroph/0403353). 

Shafieloo, A. & Souradeep, T., 2004, Phys. Rev. D, in press, (astro-ph/0312174). 

Souradeep,T. 2000, in 'The Universe', eds. Dadhich, N. & Kembhavi, A., Kluwer. 

Souradeep, T. and Hajian, A., 2004, Pramana, 62, 793. 

Spergel, D. et al., 2003, Astrophys. J. Suppl., 148, 175. 

Spergel, D., and Goldberg, D. M., (1999), Phys.Rev. D59, 103001. 

Starkman, G. Class. 1998, Quantum Grav. 15, 2529. 

Starobinsky, A. A. 1982, Phys. Lett, 117B, 175. 

Tegmark, M., de Oliveira-Costa, A. & Hamilton, A., 2004, Phys.Rev. D68 123523. 

Varshalovich, D. A., Moskalev, A. N., Khersonskii, V. K., 1988 Quantum Theory of Angular Mo- 
mentum (World Scientific). 

Vielva, P., Martinez-Gonzalez, E., Barreiro, R. B., Sanz, J. L., Cayon, L., 2003, preprint (astro- 
ph/0310273). 

Zaldarriaga, M., Seljak, U., Bertschinger, E., 1998, Astrophys. J., 494, 491. 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



- 30 - 




20 40 60 80 100 120 20 40 60 80 100 120 



l(l+l)+m+l l(l+l)+m+l 

Fig. 1. — BipoSH coefficients are linear combinations of elements of the covariance matrix. Here 
Af/} 1 (left) and Af/}' 1 (right) are plotted to show how BiPS covers the off-diagonal elements of the 
covariance matrix in harmonic space. 
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Fig. 2. — C{6) 'recovered' from a map of a statistically isotropic model matches the original C(9) 
shown as a solid line from which we made the map. The error bars shown are computed from 10 
realizations. Note that the size of error bar depends on the size of 6 bins. 
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Fig. 3. — Analytical bias for (left) a two beam window function with Wf (20,30) and (right) a 
Gaussian window function with W, (40) computed from the average C\ from 1000 realizations of a 
SI CMB map compared with ^realization^ ( the 

average K\ from 1000 realizations). This shows that 
the theoretical bias is a very good estimation of the bias for a statistically isotropic map. 




Fig. 4. — Bias corrected 'measurement' of BiPS, K£ for SI CMB maps with a best fit LCDM 
power spectrum smoothed by (left:) a two beam function W; (20, 30) and (right.) a Gaussian beam 
Wf'('iO). The cosmic error, a(ne), obtained using 1000 independent realizations of CMB (full) sky 
map matches the analytical results shown by dotted curve with triangles . This shows a much 
better fit to the theoretical cosmic variance compared to what was obtained for 100 realizations 
(Hajian & Souradeep 2003b) 
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Fig. 5. — Top: Ci of the two WMAP CMB anisotropy maps. The red, magenta and green curves 
correspond to map A, B and C, respectively. The black line is a 'best fit' WMAP theoretical 
Cg used for simulating SI maps. Blue dots are the average C\ recovered from 1000 realizations. 
Bottom: These plots show the window functions used. The dashed curves with increasing I coverage 
are 'low-pass' filter, W^{l s ), with l s = 4, 18,40, respectively. The solid lines are 'band-pass' filter 
Wf(l t ,l s ) with (l s ,l t ) = (13, 2), (30, 5), (30, 20), (45, 20), respectively. 
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Fig. 6. — Figure shows the measured values of K£ for maps A, B and C filtered with a Gaussian 
window with l s = 40. Each map appears to show the violation of statistical isotropy for some I at 
the level of ~ la. However, this is because the PDF for these models are skewed with the maximum 
probability shifted away from the mean to negative values (see Fig. 9). 
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Fig. 7. — Measured BiPS for maps A, B and C filtered with a window with l s = 30, It = 20. This 
is to check the statistical isotropy of the WMAP in the modest 20 < I < 40 in the multipole space 
where certain anomalies have been reported. ILC with a 10-degree-cut (top) has the same BiPS as 
map C (l s = 30, It = 20) which explains that the raising tail of CSSK map is because of the mask. 
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Fig. 8. — BiPS measured from 4 different WMAP maps filtered by W, (20, 45). From bottom to 
up: TOH map, ILC map, CSSK foreground free map with a galactic cut, Wiener filtered map. We 
see that Wiener filtering has the similar effect as galactic cut. (But not exactly the same) 
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Fig. 9. — Probability distribution function for K\ to k§ constructed from 1000 realizations. The left 
panel shows the PDF for the maps filtered with W z 5 (20,30) (left panel) and Wf(l s = 40) (right 
panel). The latter is more skewed, which explains the ~ la shift in K£ values shown in Fig. (6). 
The green, magenta and red (circular, pentagonal and rectangular) points represent ILC, CSSK 
and TOH maps, respectively. The smooth solid curves are Gaussian approximations. 
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Fig. 10. — The probability of the two WMAP based CMB maps being SI when filtered by 
Wf (20,30) and a Gaussian filter W, G (40). 
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Fig. 11. — Figure compares the measured values of K£ for maps A and B filtered to retain power 
only on the lowest multipoles, / = 2 and I = 3 assuming the WMAP theoretical spectrum 
WMAPbf (left) and a model spectrum that matches the suppressed power at the lowest multi- 
poles (Shafieloo & Souradeep 2004). The non zero m 'detections' assuming the WMAP theoretical 
spectrum become consistent with zero for a Cf that has power suppressed at low multipoles. 
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Fig. 12. — probability of the two WMAP based CMB maps being SI when filtered by W t G (4) to 
only keep the lowest multipoles when assuming the WMAP theoretical spectrum, Cf (left) and a 
model spectrum that matches the suppressed power at the lowest multipoles (right). 



